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A new class of analytic and parameter-free, strongly cor- 
related wave functions of simple functional form is derived 
for few electrons in two-dimensional quantum dots under 
high magnetic fields. These wave functions are constructed 
through breaking and subsequent restoration of the circular 
symmetry, and they offer a natural alternative to the Laughlin 
and composite-fermion functions. Underlying our approach is 
a collectively-rotating-electron-molecule picture. The angular 
momenta allowed by molecular symmetry correspond to the 
filling- factors' hierarchy of the fractional quantum Hall effect. 
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Two-dimensional (2D) TV-electron systems in strong 
magnetic fields have been the focus of extensive theo- 
retical investigations in the last twenty years fi]-|i~3|j. The 
reasons are twofold: (I) The early realization |l|,|2j that 
few electron systems are relevant to the physics of the 
fractional quantum Hall effect (FQHE) observed in the 
infinite 2D electron gas, and (II) The recent progress in 
nanofabrication techniques that has allowed experiments 
on 2D circular quantum dots (QD's) containing a finite 
number of electrons fJJJljl ■ 

Among the many theoretical methods for studying 
such systems, two approaches have become well estab- 
lished, i.e., exact diagonalization techniques ]i|,^p|-[lC|] 
and consideration of appropriate classes of strongly cor- 
related, analytic trial wave functions in the complex 
plane The trial wave functions proposed todate 

have been based on physical intuition, and their justi- 
fication has been inferred a posteriori through compar- 
isons with exact numerical calculations and/or with the 
phenomenology of the FQHE. 

In this paper, we use a systematic, microscopic ap- 
proach and derive a new class of strongly correlated, ana- 
lytic wave functions for the TV-electron problem in strong 
magnetic fields |HJ. Our analytic wave functions have a 
simple functional form which differs from that of the fa- 
miliar composite-fermion (CF) Q| and Jastrow-Laughlin 
(JL) (|] functions, and they are associated with a phys- 
ical picture of a collectively rotating electron molecule 
(REM). Guiding the synthesis of the states of the system, 
our approach consists of two steps: First the breaking of 
the rotational symmetry at the single-determinantal un- 
restricted Hartree-Fock (UHF) level yields states repre- 
senting electron molecules (EM's, or finite crystallites). 
Subsequently the rotation of the electron molecule is de- 



scribed through restoration of the circular symmetry via 
post Hartree-Fock methods, and in particular Projec- 
tion Techniques (PT's) p0| . Naturally, the restoration of 
symmetry goes beyond the mean-field and yields multi- 
determinantal wave functions. Earlier we demonstrated 
that this method (generalized to include in addition the 
breaking of the total-spin symmetry) can describe accu- 
rately two-electron systems in molecular |l7j and single 
Hi QD's at zero magnetic field 

In general, the symmetry-broken UHF |21 1 o rbitals 
need to be determined numerically [p~2| , p~3| l^ Jle] , |22]| . 
However, in the case of an infinite 2D electron gas 
in strong magnetic fields, it has been found that 
such UHF orbitals can be approximated by ana- 
lytic Gaussian functions centered at different positions 
Zj = Xj + iYj and forming an hexagonal Wigner crystal 
(each Gaussian representing a localized electron). The 
specific expression for these displaced Gaussians is 



i(z,Z 3 ) = (l/^) 

x exp[-|z - Zj\ 2 /2] exp[-i(xYj + yXj)} 



(1) 



where the phase factor is due to the gauge invariance. 
z = x — ly, and all lengths are in dimensionless units of 
IbV% with the magnetic length being Is = yhc/eB. 

In the case of a Coulombic finite iV-electron system, it 
has been found Jll],[l2| that the UHF orbitals arrange in 
concentric rings forming EM's (referred to also as Wigner 
molecules, WM's) The UHF results for the forma- 

tion of WM's are in agreement with the molecular struc- 
tures obtained via the conditional probability distribu- 
tions (CPD's) which can be extracted from exact numer- 
ical wave functions |) 10 2(|. For N < 4, the electrons 
are located at the apexes of a regular polygon situated 
on a single ring, while for 5 < A < 7 both the single- 
ring structure and an isomeric one with one electron at 
the center come into play. We will denote the former 
arrangement as (0, N) and the latter as (1, A" — 1). The 
electrons of the (0, N) ring are located at 



Zj = Zcxp[i27r(l - j)/N], 1 < j < N , 



(2) 



and those participating in a (1, A — 1) arrangement are 
located at 

Zi = 0; Z 3 = Zcx-p[i2n(2- j)/{N - 1)], 2 < j < N . 

(3) 

Before proceeding further, we need to expand the dis- 
placed Gaussian ([l]) over the Darwin-Fock single-particle 



1 



states. Due to the high magnetic field, only the single- 
particle states, 



(iH Hr 



N 



exp(— zz*/2) , 



(4) 



of the lowest Landau level (LLL) are needed (observe 
that the angular momentum of this state is —I due to the 
definition z = x — ly) . Then a straightforward calculation 
p7| yields 



u{z,Z)=Y t C l {Z)4> l {z) , 



(5) 



1=0 



with Ci(Z) = (Z*)'exp(-ZZ*/2)/Vn for Z ^ 0. Natu- 
rally C o (0) = 1 and C l>0 (0) = 0. 

Since electrons in strong magnetic fields are fully po- 
larized, only the space part of the many-body wave func- 
tions needs to be considered; for the symmetry-broken 
UHF determinant describing the WM, it is given by 

*tjhf = det[u(zi,Zi),u(z 2 ,Z 2 ),- ■ -,u(z N , Z N )} . (6) 

Using (^) one finds the following expansion (within a pro- 
portionality constant) 



N 

UHF 



E 



Cii (Zi)Ci 2 (Z 2 ) ■ ■ ■ Ci N (Zff) 



(i=0,...Jm=0 



N 

D(h,l 2 ,...,l N )eM-z2 ZiZ i^ ' ( ? ) 



where D(h,l 2 , ...,In) = detfz' 1 , z\ 2 , ■ 

The UHF determinant [Eq. @ or Eq. (0)] breaks the 
rotational symmetry and thus it is is not an eigenstate of 
the total angular momentum TiL = ftJ^Li U- However, 
one can restore the rotational symmetry by applying onto 
the UHF determinant the following projection operator 

M 



2ttV l 



d'y exp[ij(L — L)} 



(8) 



where TiL = h^2 i=1 li are the eigenvalues of the total 
angular momentum. 

It is advantageous to operate with Vl on expression 
(0) , which is an expansion in a basis consisting of prod- 
ucts of single-particle eigenstates with good angular mo- 
menta U. Indeed in this case the projection operator acts 
as a Kronecker delta: from the unrestricted sum (Q), it 
picks up only those terms having a given total angular 
momentum L. As a result, after taking into consider- 
ation the specific electron locations (||) associated with 
the (0, N) WM, one derives (28) the following symmetry- 
preserving, many-body correlated wave functions (within 
a proportionality constant), 



e nw 

o<h<h<—<iN y 



' 1 I~T s * n 

,l<i<j<N 



N 



(k - l 3 ) 



N 



D(k,l 2 , -Jn) exp(- ^2 z i z i/ 2 ) 



(9) 



In deriving (^) , we took into account that for each deter- 
minant D(li,l 2 , ■■■,In) in the unrestricted expansion (|7]) 
there are N\ — 1 other determinants generated from it 
through a permutation of the indices l 2 , In}] these 
determinants are equal to the original one or differ from 
it by a sign only. In the case of an (1, N — 1) WM, the 
corresponding correlated wave functions are given by, 



l 2 + ... + l N= L f N \ 

E 

0<1 2 <1 3 < — <1 N \i=2 / 



\ ! J~[ sin 

\2<i<j<N 



N-l 



x D(0,l 2 , ...,l N )exp(- 



N 

E 



Zl z*/2) 



(10) 



We call the correlated wave functions [Eq. (^) and Eq. 
(|l0|)1 the electron-molecule wave functions (EMWF's). 
We stress that the EMWF's have good total angular 
momenta, unlike the UHF determinant from which they 
were projected out. The projection operator (pft acts on 
a single UHF determinant, but yields a whole rotational 
band of the WM. The states in this band are those with 
the lowest energy for a given angular momentum L, and 
in addition they are purely rotational, i.e., they carry no 
other internal excitations; in analogy with the custom- 
ary terminology from the spectroscopy of rotating nuclei 
p6| , p9[ , we designate this band as the "yrast band" . 

Furthermore, if instead of electrons the displaced 
Gaussians ([l]) describe bosonic particles forming a 
molecule, the corresponding ]lq ] many-body correlated 
wave functions will be given by expressions similar to 
Eq. (^) and Eq. (|To|), but with the following two impor- 
tant differences: (I) The product of sine functions will 
be replaced by a sum over cosines, and (II) The determi- 
nants D(l\, l 2 , In) will be replaced by permanents |3C| l 
P(li,l 2) ...,l N ) = pcrm^ 1 , z l 2 2 , z l £]. 

Among the properties of the EMWF's specified by Eq. 
@ and Eq. ([!(]), we mention the following: 

1) The EMWF's lie entirely within the Hilbert sub- 
space spanned by the lowest Landau level and, after ex- 
panding the determinants |^8| , they can be written in the 
form (within a proportionality constant), 
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N 

*£[*]=i^[aM-£*i*i/2), (11) 
i=l 

where the [z] 's are order- L homogeneous polynomials 
of the Zi's. 

2) The polynomials P^\z\ are divisible by 

Pv\A= II ( 12 ) 

l<i<j<N 

namely P^ [z] = P^ [z] [z] . This is a consequence of 
the antisymmetry of [z] . Py [z] is the Vandermonde 
determinant D(0, 1, AT). For the case of the lowest al- 
lowed angular momentum Lq = N(N — l)/2 (see below), 
one has P^ [z) = Py[z\, a property that is shared with 
the Jastrow-Laughlin [pj and composite-fermion [gj trial 
wave functions. 

3) Upon the introduction of the Jacobi coordinates, 
the center-of-mass separates from the internal variables 
in complete analogy with the exact solution. 

4) The coefficients of the determinants [i.e., products 
of sine functions, see Eq. (|| and Eq. flTc|)] dictate that 
the EMWF's are nonzero only for special values of the 
total angular momentum L given by, 

L = N(N — l)/2 + Nk, k = 0,1,2,3, ... , (13) 

for the (0, N) configuration, and 

L = N(N — l)/2 + (JV— l)fc, fc = 0,l,2,3,... , (14) 

for the (1,N — 1) one. The minimum angular momen- 
tum Lq = N(N — l)/2 is determined by the fact that the 
D determinants [see Eq. (^) and Eq. (|l0|)] vanish if any 
two of the single-particle angular momenta U and lj are 
equal. In plots of the energy vs. the angular momenta, 
derived from exact-diagonalization studies ]5|-p0[|, it has 
been found that the special L values given by Eq. fll3] ) 
and Eq. ( |l4| ) exhibit prominent cusps reflecting enhanced 
stability; as a result these L values are often referred to 
as "magic angular momenta". We stress that the an- 
gular momenta associated with the EMWF's correspond 
precisely to the magic L's of the exact-diagonalization 
studies |nj. In the thermodynamic limit ^,|), one can 
relate the total L to a fractional filling through the rela- 
tion v = N(N - 1)/(2L), and thus the EMWF angular 
momenta ( |l3| ) and ([L4|) correspond to all the fractional 
filling factors associated with the FQHE, including the 
even- denominator ones, i.e., v = 1, 3/5, 3/7, 5/7, 2/3, 
1/2, 1/3, etc... 

5) For the case of two electrons (N — 2), the EMWF's 
reduce to the Jastrow-Laughlin form, namely 

pim= n (*-*i) L > ( i5 ) 

l<i<j<N 

where L = 1, 3, 5, ... However, this is the only case for 
which there is coincidence between the EMWF's and the 



TABLE I. The Qg[z] polynomial associated with the 
EMWF's and the JL functions (The [z] polynomials are 
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JL wave functions. For higher numbers of electrons, N, 
the EMWF polynomials P^[z] (apart from the lowest- 
order Vandermonde P^ [z] ones) are quite different from 
the corresponding JL or composite-fermion polynomials. 
In particular, the familiar factor Yli<i<j<N( Zi ~ z i) 2% '> 
with p an integer [^Q, (which reflecFs multiple zeroes) 
does not appear in the EMWF's (see, e.g., Table I which 
contrasts the Qg[z] polynomials corresponding to the 
EMWF's and JL functions). 

6) For the case of three electrons (N = 3), after trans- 
forming to the Jacobi coordinates z = (zi + Z2 + zq)/3, 
z a = (2/3) 1 / 2 ((z 1 + z 2 )/2 - z 3 ), z fc = (zx - z 2 )/V2 
(and dropping the center-of-mass exponential factor), the 
EMWF's can be written as (again within a proportion- 
ality constant), 

<f> 3 L [z a , z b ] = [(z a + iz b ) L - (z Q - iz b ) L ] 

x exp[(-l/2) (z Q z* + z h zl)] , (16) 

with L = 3to, m = 1, 2, 3, 4, ... being the total an- 
gular momentum. Again the wave functions $^ [z a , z b ] 
are very different from the three-electron JL ones; e.g., 
they are nonvanishing for even m values, unlike the three- 
electron JL functions. However, the ^\[z ai z b \s coincide 
with the functions \m, 0) derived in Ref. jl|. We notice 
that, although it was found Jl],|32| that these wave func- 
tions exhibited behavior expected of fractional quantum 
Hall ground states, the generalization of them to a higher 
number of electrons did not follow. 

Several publications |^|||[33| have applied the compos- 
ite fermion picture (the JL functions are a special case of 
the CF's) to single QD's in strong magnetic fields. In par- 
ticular, it has been shown Q that CF wave functions can 
be constructed with angular momenta coinciding with 
the magic ones. However, it has also been found S that 
several discrepancies exist, i.e., some of the larger magic 
angular momenta are not reproduced by the CF picture. 
As a consequence of the above, the REM description with 
the EMWF's derived here offers a natural alternative for 
interpreting the physics of electrons in QD's in high mag- 
netic fields. This proposition is further supported by in- 
spection of the overlaps between the EMWF's and the 
exact many-body eigenstates, and their comparison with 
the corresponding overlaps for the JL states; see Table 
II, where in some instances (i.e., N — 4, L — 10 and 14) 
we list energies of the EM, CF, and exact states instead 
of the overlaps. Indeed the agreement between the EM 
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TABLE II. Overlaps, >/({^ |^> (V# |v^» 1/2 , of 

EMWF's (0's) and JL functions (</>'s) with the correspond- 
ing exact eigenstates (tp's) for various values of the angular 
momenta L. Recall that the angular momenta for the JL 
functions are Ljl = N(N — l)m/2, with m > being an odd 
integer. Bottom: Energies of EMWF's compared to CF and 
exact-diagonalization results. Energies in units of e 2 /kIb, (k 
is the dielectric constant). 



OVERLAPS 


L 


EMWF 


JL 




N=3 


9 


0.98347 


0.99946 a 






15 


0.99473 


0.99468 a 






21 


0.99674 


0.99476 a 






27 


0.99758 


0.99573 a 






33 


0.99807 


0.99652 a 






39 


0.99839 


0.99708 a 




N=4 


18 


0.92937 


0.97880 






30 


0.96742 


0.94749 






42 


0.97366 


0.95561 






54 


0.97623 


0.96815 




ENERGIES 


L 


EMWF 


CF 


EXACT 


N=4 


10 


1.78510 


1.78537 b 


1.78509 




14 


1.50955 


1.50222 b 


1.50066 



a From Ref. §. 

b From Table V of Ref. Eil 



states and the exact ones is of comparable quality as in 
the case of the CF and JL wave functions. 

In summary, we have developed a new class of analytic 
and parameter- free, strongly correlated wave functions 
of simple functional form, which accurately describe the 
physics of electrons in QD's under high magnetic fields. 
The thematic basis of our approach is built upon the in- 
tuitive picture of collectively rotating electron molecules, 
and the synthesis of the many-body EMWF's involves 
breaking of the circular symmetry at the UHF level with 
subsequent restoration of this symmetry via a projection 
technique. While we focus here on the strong magnetic- 
field regime, we note that the REM picture unifies the 
treatment of strongly correlated states of electrons in 
QD's over the whole magnetic-field range p"2| , p2| , p6| . We 
also remark that our analysis, aimed here mainly at treat- 
ing finite electron systems (i.e., QD's) with an arbitrary 
number of electrons, points to the remarkable conclusion 
that the observed FQHE hierarchy of filling factors may 
be viewed as an experimental signature of the yrast band 
(see above) of the REM. 

This research is supported by the U.S. D.O.E. (Grant 
No. FG05-86ER-45234). 
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